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Abstract 

The purpose of this paper is to come up with a framework that "converts" existing 
concepts from configuration space to ordinary one. This is done by modehng our 
universe as a big "computer" that simulates configuration space. If that "computer" 
exists in ordinary space and is ran by "classical" laws, our theory becomes "classical" 
by default. We have first applied this concept to a version of quantum field theory in 
which elementary particles have size (that is, a theory that does not yet exists). After 
that, we have also done the same with Pilot Wave model of discrete jumps, due to 
Diirr et el. 

1. Introduction 

When we think of non-relativistic quantum mechanics of a single particle, we can re-interpret 
it as a classical wave mechanics: a field ip "evolves" according to "classical" Schrodinger's 
equation. We do not have to think of it as probability. It simply happens that the probability 
of the "collapse" of that "classical" field is proportional to its amplitude square. Furthermore, 
Bohm's Pilot Wave model can provide a mechanism by which the probability of finding a 
particle at any given place is proportional to I'i/'p (even though in case of single particle there 
is no "collapse" since there is no entropy required for decoherence). 

However, when we introduce more than one particle, this qualitative picture completely 
changes. The function is no longer a function on a position space but rather on a config- 
uration space, as it assigns complex amplitudes to different configurations of particles; and, 
in quantum field theory it becomes a function on Fock space, where the numbers of particles 
are no longer fixed. Since we are used to think of a "field" as a function strictly on M^, 
we can no longer view probability amplitudes as one. This forces us to take the word of 
"probability" in the term "probability amplitude" more seriously, which brings us back to 
the paradox of its complex value. 

The ultimate answer to this question is to restore the usual three dimensional space. 
It is possible to do so by the following argument. Suppose we were living in a classical. 
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three dimensional world. In that world, we could have designed a computer program that 
simulates configuration space, probability amplitudes, and everything else we have to deal 
with in quantum mechanics. Now, there is no such thing as " computer program" . In reality, 
computer is made out of particles, and what it shows in the screen is a result of complex 
interaction of particles. This implies that there are three dimensional processes taking place 
that "simulate" non-existant, multidimensional configuration space. In this work we will 
come up with one such three dimensional classical device that does it. 

We then take advantage of this philosophy and notice that we can "write a computer 
program" not only for computing amplitudes of quantum states, but we can also write a 
computer program for Pilot Wave model as well. For the purposes of this paper, we will 
"write a program" for the model with Stochastic jumps proposed by Diirr et el (see [1]). We 
chose this specific approach because it invokes discrete quantum states, and, as we shall see, 
discrete states are easiest to model that way. But, at the same time, our "computer" is ran 
by classical laws that are smooth and deterministic. This cures the model from the violation 
of determinism due to stochastic nature of originally proposed jumps. In principle, the same 
can be done for any other Pilot Wave model. But, since the latter is a lot more complicated, 
we leave other Pilot Wave models for future research. 

It has to be pointed out, however, that while our "computer" is ran by "local" signals, 
their speed of propagation is much faster than the speed of light. After all, it is not possible 
to build a computer that simulates signals faster than the ones in the world in which it 
lives. After all, any signal that computer simulates is a result of a set of signals between 
its particles. However, since our intuition does not demand relativistic covariance, I don't 
regard it as a big problem. The only things our intuition does demand are (non-relativistic) 
locality and determinism, and the mechanism by which our "computer" operates is both 
local and deterministic. We then argue that the appearance of relativity in the lab is only a 
result of the specifics of our Hamiltonian. 

On the first glance it might sound like cheating. After all, we could use this argument 
for anything and everything. Suppose, for example, we didn't like Coulumb's law, and we 
liked the l/r^° a lot better. We could then model computer, based on l/r^° that simulates 
on its screen. To our defense, we will point out that we always arrive at the same 
problem when we are trying to do something we don't know. Suppose, for example, Newton 
didn't have the numerical information that he used in analysis of Kepler's laws. Then, if he 
would come up with the same law of gravity, it would look very arbitrary. 

From the latter point of view, the ultimate reason why our theory looks arbitrary is 
simply that we don't have any information of what happens on the small scale, which is 
what forces us to be creative. While this is not a pleasant thing to hear, it is a lot better 
than saying that our very classical intuition on that scale is wrong. Of course, that might 
be the case, but it doesn't have to be. The purpose of this work is to show one way in which 
the processes on small scale might look classical. This, of course, is just one way out of many 
other alternatives, and I do not claim this to be the truth. 
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2. How wave function is "encoded" 

According to our model, the configuration space (or Fock space in case of quantum field 
theory) is discretized. Every discrete state is represented by a subset of R^. More specifically, 
we have a fixed configuration of particles in R^. By specifying a specific subset of R^ we 
are choosing to "look at" the particles that are inside that subset, while "ignoring" all other 
particles. This specifies the quantum state. Now, if we are going to alter the subsets of 
we are "looking" at, this will result in a perception that the particles are either created or 
destroyed. Now, if we make sure that every such subset of R^ is "dense enough" (that is, 
if small enough neighborhood of every point intersects all of the subsets in our list), then 
this will allow for a perception of continuous processes, such as motion of particles, on small 
enough scales. 

To see how this model works, consider a simple example. Suppose the configuration 
space has only 3 elements. We assume that we have one single configuration of particles, 
none of which can move (and, therefore, their common reference frame is a "preferred one"). 
So, for example, we have 5 electrons, 4 protons, and 8 photons. We assume a toy model 
of only having one space and one time dimension. So, the x coordinates of 5 electrons are 
4.89, 1.74. 6.95, 5.26 and 8.31. The x coordinates of 4 protons are 3.76, 2.38, 6.11 and 4.75. 
Finally the x coordinates of the photons are 2.98, 4.64, 7.22, 7.23, 8.11, 4.87, 1.39 and 3.68. 

Now we will break set R into three subsets, 5*1, 5*2 and 5*3. The set Si consists of all 
numbers whose second digid after the "dot" is neither divisible by 2 nor by 3. The set 5*2 
consists of all numbers whose second digit after the "dot" is divisible by 2 but not by 3. 
Finally, the set S2, consists of all numbers whose second digit after the "dot" is divisible by 

3. In general, the number of such sets is the same as the number of points in configuration 
space (which, of course, is a very large number). As stated previously, in order to allow the 
combination of creation and annihilation to create an illusion of continuous motion, these 
sets have to be defined in such a way that they intersect small enough neighborhoods of each 
point, and the measures of these intersections are roughly the same. Apart from that, the 
details of their definitions are not very important. 

Now, to 5*1, 5*2 and 5*3, we associate quantum states |si >, \s2 > and |s3 >, respectively. 
These states consist of configurations of particles whose coordinates happen to be elements 
of Si, S2 and 5*3, respectively. Thus, \si > consists of two electrons, at locations 6.95 and 
8.31, two protons, at locations 6.11 and 4.75, and two photons, at locations 8.11 and 4.87. 
The \s2 > consists of one electron at a location 1.74 (5.26 is disqualified since 6 is divisible 
by 3), one proton at a location 2.38, and four photons at locations 2.98, 4.64, 7.22 and 3.68. 
Finally, |s3 > consists of 2 electrons at locations 4.89 and 5.26, 1 proton at a location 3.76, 
and two photons at locations 7.23 and 1.39. 

Now, we introduce a complex valued field ip on our space R. We come up with a local 
dynamics of ip in such a way that at the equilibrium ipix) ~ ip{x') as long as both x and 
x' are elements the same set Sk, for some k. If such is the case, then we can simply define 
the amplitude of |sfc > as a value of ip{x) for any x e Sk- At the same time, in order for 
the quantum amplitudes to evolve, that probability amplitude changes with time. This, 
of course, violates the relativistic covariance, and it logically corresponds to the violation 
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of relativistic covariance in the "standard" definition of configuration space. However, this 
new picture allows us to continue to view quantum mechanics as a classical field theory in 
ordinary space, albeit non-relativistic one. 

It should be pointed out that an element of the set Sk that "represents" \sk > does not 
have to coincide with the locations of any of the particles. Thus, ip{57.11) is approximately 
equal to the amplitude of |si >, even though the point x — 57.11 is more than 47 units away 
from any of the particles in configuration |si >. On the other hand, ■0(6.94) approximates 
the amplitude of \s2 > and not \si >, despite the fact that x = 6.94 is only 0.01 units away 
from the electron at x = 6.95 which is part of |si > state. In other words, despite the fact 
that the dynamics of ijj is local, at the equilibrium ijj might change very rapidly in space, 
and yet have similar values at the points that fall into the same state, no matter how far 
away they are. 

The particles, themselves, do not move. But whenever a "measurement" is performed, 
-0 collapses into an "extreme state" where it is 1 at all x E Sk and everywhere else. Thus, 
all the particles whose coordinates are not elements of Sk are "hidden from view". Thus, 
if multiple measurements are performed, different particles get "hidden" each time. If we 
have two electrons very close to each other, and first electron was hidden at ti but not at t2 
while the second one was not hidden at ti but became hidden at t2, that might lead to an 
appearance that electron have moved from one location to another, provided that t2 — ti is 
small. In reality, however, both of these electrons were always stationary in their respective 
locations, and always existed as separate particles. 

Let us illustrate a point with our example, of a configiuation space consisting of three 
points. Suppose we performed two separate measurements, at t = ti and t = t2. We have 
found that, aX t = ti, ip{37.91) = 0(69.32) = and -0(107.13) = 1; on the other hand, at 
t — t2 we got 0(25.61) = 0(273.83) = and 0^(415.72) = 1. Prom this we know that at time 
t — ti our system had collapsed to |s3 > and at the time t — t2 it had collapsed to |s2 >■ 
Thus, at ti we had one electron at a location 1.74, one proton at a location 2.38, and four 
photons at locations 2.98, 4.64, 7.22 and 3.68. Then, at ^2 we had 2 electrons at locations 
4.89 and 5.26, 1 proton at a location 3.76, and two photons at locations 7.23 and 1.39. 
As far as most particles are concerned, we wrongly interpret it as them being created and 
annihilated. Furthermore, we also wrongly decided that one of the photons have "moved" 
from location 7.22 to 7.23. In reality, none of the particles were created or destroyed, and 
none were "moving" either. Both illusions were accomplished through particles hiding from 
our view. 

It should also be pointed out that, no matter how we define different Sk, and no matter 
what the stationary particle configuration is, we have a "global" law that if at least one of the 
particles "moves", then all have to "move". That is due to the fact in order for at least one 
particle to move, we have to switch states. But, since every point in space is "reserved" for 
only one state, once the states were switched, none of the particles can be found at any point 
"occupied" by the old state. In order for this not to contradict our experience, we have to 
make sure that any Sk intersects with small enough neighborhood of any given point. In this 
case, up to some coarse-graining, every point will be approximately represented in any given 
state. Furthermore, we would like to make sure that the intersections with neighborhoods 
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wc have just mentioned have roughly the same measure. This allows us to avoid unwanted 
variations of "densities" of the particles on larger scales. 

It should also be noticed that, due to the fact that particles do not move, in our example 
we can not obtain any states besides the three that we listed. The reason in reality we have 
every possible configuration of particles is simply because we have a lot more particles, and a 
lot more S-s. Once these numbers are large enough, we can statistically expect to generate an 
approximation to any conceivable state, in at least one of these S-s. Now, statistically we can 
also expect to generate several S-s that simulate the same state. But, again from statistical 
point of view, we expect this overcounting to be approximately the same for different S-s, 
which allows us to expect the probability amplitudes to be consistent with quantum field 
theory. 

The main challenge of this work is to obtain the global picture we have just outlined 
thorough the local dynamics of ip and some other fields. In particular, we would like our 
local theory to predict that, at any given moment in time, ip can vary a lot within any 
small neighborhood, and yet it is constant within each set, despite the fact that our sets are 
"spread out" throughout the whole universe. We accomplish this goal by introducing two 
stationary fields, Ui{x) and U2{x) and comming up with a dynamics that assures us that each 
point in space receives signals with frequency Ui{x) and emits the ones with frequency U2{x). 
In this case, the points with similar values of ui are "listening" to the same wavelength. 

To illustrate the way it works, suppose there are three people. A, B, and C. People A 
and B are neighbors, while a person C lives at the opposite part of the world. A person 
A and a person C are both listening to the radio station 1, while a person B is listening 
to radio station 2. If these people are completely isolated from any sources of information 
other than their respective radio stations, then the views of a person A will be the same as 
the views of the person C and not person B, despite the fact that person B is very close to 
a person A, while person C is far away. At the same time, the picture is completely local, 
since the mechanism by which radio waves spread is a local one. 

Of course, in case of the radio we have amplifiers that assure us that signals don't get 
weakened in space. Since in our case we don't have these, we instead assume that the 
universe is compact (for example, a large sphere). In light of the fact that we admit that 
we violate relativity, we also assume that the speed of radio waves is much faster than the 
speed of light; in fact, it is so large that they can circle the universe within a very short time. 
In light of this, our waves circling the universe several times and, as a result, it is no longer 
important where they were emitted; their intensity, on average, is the same everywhere. Due 
to the compactness of the universe, their intensity does not go to zero. 

Finally it is improtant to address the unnatural behavior of Ui and U2. First of all, as 
mentioned earlier, they have to be very fast varying. Secondly, in order for the technicalities 
of the theory to work, they have to be nearly integer most of the time. Both of these 
assumptions are very unnatural. In order to make the theory more natural, we postulate 
vi and V2 as basic fields, both of which are slowly- varying, differentiable functions (and, of 
course, non-integer most of the time). Wc then identify Ui and U2 with g{vi) and g{v2), 
respectively, where g{x) has the properties that we would like ui and U2 to posess. For 
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example, we can set 



g{x) — — tan ^{l{m sin x — [m sin x\)) (1) 

k=—n 

where /, m and n are very large integers. The fact that / is large assures us that each term in 
the sum is either close to 7i/2 or — 7r/2. Thus, the factor 2/71 assures us that each term is close 
to either —1 or 1. Thus, the overall sum is close to any integer between —2n — 1 and 2n + l. 
Thus, in order to get the desired number of points in configuration space, n has to be equal 
to that number divided by 2. Finally, the fact that we have used m sin x—\m sin x\ instead 
of x itself assures us that even when Vi and V2 are smoothly increasing (or decreasing), the 
functions Vi and V2 still quickly fluctuate between all integer values between —2n and 2n, 
as desired. The fact that m is large implies that these fluctuations are quick. As a result, by 
setting ui — g{v\) and U2 — (7(^2) we obtain desired behavior for ui and U2 while postulating 
vi and V2 to be of a lot more natural form. 



3. The Hamiltonian 

In the previous section we have explained the way we define the complex amplitudes in 
a configuration space as a function in an ordinary, three dimensional one. We have also 
defined the concept of evolution of our states, and have shown how apparent creation and 
annihilation of the particles is possible despite the fact that in reality their number is fixed. 
We have likewise shown how appearance of their motion is possible while they are stationary. 
All of these are kinematical concepts. In this section we will move to dynamics. In other 
words, we will explain what we mean by the Hamiltonian that guides the evolution of a wave 
function. 

Again, we are working with a toy model where configuration space has only three points. 
But this time, in addition to sets 6*1, S2 and ^3, we will introduce the sets Ti, T2 and T3. 
The sets Si, S2 and are defined as before; on the other hand, Ti, T2 and T3 are defined 
based on the first digit after the dot rather than the second one; the rest of the definition 
is the same. So, for example, 3.14 is an element of yS'2 fl Ti. For any x & Sk C\ Ti, we 
would like e*^^^^ to represent the probability amplitude of the transition from \sk > to \si >. 
Thus, e^^^^'^^^ is a probability amplitude of transition from \s2 > to |si >, while e''^^^'^^^ 
is a probability amplitude of transition in the opposite direction. The unitarity condition 
implies that H(3A1) = -H{3.U). 

Now suppose we only have three electrons, and they are located at 5.31, 5.32, and 7.83. 
Thus, the state |si > corresponds to one electron at 5.31, the state |s2 > corresponds to one 
electron at 5.32 and the state |s3 > corresponds to one electron at 7.83. From completely 
local point of view, none of the states can transition to each other, since none of the particles 
are "touching" . Given that, from our setup, it is always bound to be the case, we instead say 
that the transition can occur if two particles are close enough (we will discuss this issue in 
more detail shortly). Thus, we would like the probability amplitudes of transitions between 
|si > and |s2 > to be large, while the probability amplitudes between |s3 > and any of the 
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other two states to be small. This means that we would like e^^i"^-^^) and e'-^^^'^^^ to be large, 
while having e^^^^-^^\ ^iH{2m) ^ ^iH{5.i3) ^iH{9.35) ^^^n 

Now, suppose we add another electron at 8.72. Then the probability amplitude of the 
transition between 5*1 and 5*2 will become small, since there is no electron near 8.72 in 5*1. 
This phenomenon can be explained by expressing H as an integral of a local function h{x), 

J h{x)Sx (2) 

In the above expression, e^^^^^ represents a probabilityamplitude of an imaginary quasi-local 
transition, with some specified "small" scale (in our case, lets say that scale is 0.2) and e*^ 
on the other hand represents the probability amplitude of the transition between two global 
states. Since exponent of the sum is a product of exponents, simple algebra tells us that e*^ 
is a product of e^^, as desired. 

In our case, despite the fact that the probability amplitude of transition between 5*2 and 
(S'a are small, e*'^*^^'^^-' is large. After all, as far as neighborhood of 2.43 is concerned, there 
are no particles in either state and the probability amplitude of transition between empty 
states is large. At the same time, however, we assign to H pointwise values too, not just to 
h (namely, for any x e Sk, H{x) ^ H{S)). In light of this, e"-"^"^-^^^ is small, despite e''^^^-^^^ 
being large. This is due to the fact that e''*^^-^^) and e*'**^^'^^) are small (due to the electrons 
at these points). Both contribute to e*^^^'^'^^ but not to e*'*^^'^^\ This is despite the fact 
that both H and h are defined at a point (namely, at 2.43). In other words, as far as h is 
concerned any given point is "listening" to (a part of) its neighborhood, and as far as H is 
concerned it is "listening" to the points far away, as well. 

The "local" mechanism of performing the above integral is modeled as follows: Every 
point emits a signal with a certain frequency based on the values of fields ui and U2. We 
recall from previous section that Ui and U2 are close to integer most of the time, with an 
exception of transition regions that are necessary to make sure they are continuous. Their 
integer approximations correspond to the points in a discrete configuration space (thus, 
Ml (5. 12) = 2 and ti2(5.12) = 1). Now the frequency of the emitted signals is a function of 
both ui and U2 (we will call it uj{ui,U2)). In light of discreteness, we can define ui in such a 
way that a;(iti, U2) ~ u;(w'i, u'^ if and only if hath ui ~ u\ and U2 ~ u'2 hold. 

Now, the signals that are to be received with these resonance frequencies communicate 
the local value of h. The only other points that can receive that signal are the ones that 
have the same ui and U2. In other words, they are in the same Si n Tj as the original point. 
Thus, a signal emitted at 6.37 can not be received at points such as 4.97, 6.38 or 9.73. In 
case of our example, the only points that can receive that signal are the ones that fall in 
the same category based on both digits. The examples of the latter are 5.37, 9.91, 6.95, etc. 
Once these signals are received, they are "converted" into non-oscillating values (based on 
the mechanism discussed in chapters 5 and 6). As a result the value of H is determined by 
a sum of the amplitudes of all of the signals that any given point receives. Thus, if we allow 
the signals propagate much faster than the speed of light, and claim that they can circle a 
universe within a very short time, then H will approximate the integral of h. The details of 
this will be done in Section 6. 
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Another thing that needs to be mentioned is that, despite quasi-locahty, h is uhimately 
expressed locally due to our clever way of making particles "look like" they have non-zero 
size. More precisely, when we say that a particle has a size, we mean that there is some finite 
differentiable function of space, which represents the density of matter, which happened to 
be a differentiable approximation to 5-function. Now, in order to accommodate our desired 
picture, we will continue to view the particles themselves as points, but this density function 
will be a short-range field that our point particles produce. Each kind of particle has its 
own density field: the density field of electron is p^, the density field of a proton is Pp, the 
density field of positron is Pe+, etc. The quasi-local interaction between the particles is due 
to the local interaction of their density fields. 

There is one subtlety here. Namely, despite the fact that h "looks" at small but finite 
neighborhoods, we do not want it to look continuous on that scale either. For example, if 
we consider distances of the order of 0.2 to be "small" then e^K^-^"^) is large, while e*'^^^-^^) is 
small. After all, the "local" electron at 5.32 is only relevant for the former but not for the 
latter. In order to account for this, we have to make sure that p has similar behavior. So, for 
example, pe(7.82) = 0, despite the electron at 7.83. At the same time, pe(7.79) is non-zero 
due to that same electron. The way the above is accomplished is that each point is "tuned 
on" to its own frequency. Thus, when a point particle occupies a certain location, it sends 
signals to its neighboring points "telling" them to increase p. But only the points that are 
"tuned in" to the same frequency are able to receive these signals. Thus, they are the only 
ones that actually increase p. 

Now, in order for the interaction to occur, there has to be two kinds of p-s: pi, which 
corresponds to the particles that went into the process and p2, which corresponds to the 
particles that are coming out of the process. Every point has receptors for both of these 
fields; the latter are set up based on u\ and U2 fields, as usual. The field pi is being received 
based on ui which, in our case, corresponds to the second digit after the dot; on the other 
hand, p2 field is being received based on U2 which, in our case, corresponds to the first digit 
after the dot. Thus, at a point 7.95, p2 is large, while pi is zero. On the other hand, at 7.79, 
Pi is large while p2 is zero. 

In order for e^^ to be large, pi and p2 have to correlate. For example, pi(5.12) is large 
due to the electron at 5.32, and p2(5.12) is large due to the electron at 5.31. Thus, e*'^^^-^^^ 
is large. On the other hand, pi(5.14) is small since electron at 5.31 no longer counts. On 
the other hand, p2(5.34) is still large due to an electron at 5.31. Thus, e^^^^'^^^ is small. As 
one can see, this reproduces the behavior of e"^ we have outlined earlier; but this time the 
mechanism is completely local. 

In case of different kinds of particles, e''' depends on correlation of these densities. It 
closely correlates to Feynman diagrams. For example, since we know that there is a process 
e — >■ Wv, we also know that if pie, P2W ^-nd p^u are all large, there has to be a possibility 
for e*'^ to be large as well. The reason I say the word "possibility" is that if there arc other 
particles they might prevent this from being the case. For example, suppose that on top of 
the above, we also have a proton field p^; p2p happened to be large while pip is small. Then, 
the local transition would imply a production of proton from weak interaction. Since we do 
not want the latter, we set e*'* to be small despite the fact that it would have been large 
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judging from Pv and pw alone. 

It is easy to see that the above imphes that the dependence of h on p-fields serves as a 
replacement of vortexes in Feynman diagrams. This, of course, requires rewriting quantum 
field theory in terms of particles that possess size. A lot of it is very non-trivial. For example, 
Feynman diagrams are generally solved in momentum space; we will now have to redo it for 
the position space. This might be less obvious than simple Fourier transform since some of 
the Feynman rules, such as momentum conservation at the vortexes, are designed specifically 
for momentum space. Furthermore, in order to allow for propagation of particles, we would 
have to look at the derivatives of the densities. These derivatives will tell us in what direction 
the particles are displaced from each other, and thus enable us to include that information 
into a kinetic part of the Hamiltonian. 

Of course, both of the above points call for some machinery that is very different from 
what is normally being done in quantum field theory. It requires a proposal of a specific func- 
tion h{ps, pp, ...), followed by a rigorous mathematical proof that this specific function will, 
indeed, approximate standard quantum field theory (both kinetic and interaction terms). 
This is a very difficult project and is far beyond the scope of this paper. For purposes of this 
paper, we will pretend that this was already done, and show how determinism and locality 
is restored once this is the case. 

4. Desired global picture at the equilibrium 

In the past three sections we have outlined qualitatively what we are trying to accomplish. 
Now it is time to turn it into more quantitative work. In the next three sections we will 
quantify the above by means of the local theory of signal propagation. However, since the 
picture of configuration space is a global one, we would like to go from top to bottom. In 
this section we will describe the global picture that we would like to obtain at the end of the 
day. Then, in the subsequent sections, we would attempt to come up with a "local" theory 
that approximates this global picture at the equilibrium. 

Since we would like to use emission and absorption of waves to arrive at the picture, we 
know that anything that is subject to oscillation changes both magnitude as well as a sign 
in time. Thus, the only physically meaningful interpretation of i/j is the amplitude of the 
oscillation. However, the amplitude is always positive and real. In order for ip to have both 
imaginary and negative values, we instead define four distinct fields ipi, ip^i, ipi and ip^i and 
define ip as 

ijj ^ijji- + iijji - (3) 

Each of these four fields has its own oscillation, and each corresponds to the amplitude of 
that oscillation. Thus, each of these four fields is positive and real, but the coefficients of —1 
and ±i allow ip to have all possible complex values. One notices that there are many choices 
of the values of these four fields that give the same value of ip. In particular, ip is covariant 
under the transformations 

Ipl^ipl + Xr] 1p2^1p2- Xr (4) 
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and 

A^A + Xc] i>-i - Xc (5) 

For reasons that will become clear later, we will use this freedom and claim that all four 
fields grow in time, and their growth cancels each other in such a way that the resulting field 
%Ij behaves in expected fashion. 

Let us now come up with a dynamics of these fields. We will begin by describing 
dynamics in terms of regular, non-local, configuration space (or, in case of quantum field 
theory, Fock space), which we will "convert" into our usual space shortly thereafter. To 
distinguish elements of configuration space (or Fock space) from the ones of usual space, 
we will denote the former by S and the latter by x (notice that only the elements of the 
regular space have a vector sign). We will assume that the configuration space is discrete. 
Therefore, the most general dynamics for ip ui a. configuration space is given by 

m^^Y.^„(s,'S,)M.,) (6) 

I 

Just like we did with ■0, we will break H into four pieces: 

H^Hi-H_i + iHi - iH_i (7) 
By substituting the values of ip and H, the right hand side of the evolution equation becomes 

J2 iH{Sk, Si)^{Si) - J2^H^{Sk, Si)MSi) + H.,{Sk, Si)^.i{Si) 
I I 

+ H,{Sk, Si)ip.i{Si) + H_,{Sk, Si)'il^i{Si) - H^{Sk, Si)ip.,{Si)- 
- H_,{Sk, Si)MSi) - H,{Sk, Si)^,{Si) - H_,{Sk, Si)i^^i{Si)+ (8) 
+ iH,{Sj,, Si)^i{Si) + iH_^{Sk, Si)^^i{Si) + iH,{Sk, Si)^,{Si)+ 
+ iH_i{Sk, Si)ij_i{Si) - iHi{Sk, Si)^_i{Si) - iH_i{Sk, Si)'4)i{Si)- 
- iHi{Sk, Si)p_i{Si) - iH_i{Sk, Si)ilJi{Si) 

By equating it to dip{t, Sk)/dt, and remembering that each of the i/ji, ij^i and ijj^i 
has to be positive, we read off the dynamic equations for these four fields: 

dtMSk,t) = J2iHi{Sk,Si)^..{Si,t) + H_,{Sk,Si)ij,{Si,t)+ 
I 

+ H,{Sk, Si)ij_i{Si, t)) + H_,{Sk, Si)MSi, t) 

dtiP-i{Sk, t) = ^(i/i(^fe, Si)ilJi{Si, t) + H_,{Sk, Si)ij.i{Si, t)+ 
I 

+ H,{Sk, Si)MSi, t) + H^,{Sk, Si)i/j-i{Si, t)) (9) 

dMSk,t) = Y,iHiiSk,Si)MSi,t) + H_^{Sk,Sl),p_^{Sl,t)+ 
I 
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+ HiiSk, Si)i;_i{Si, t)) + H_i{Sk, Si)A{Si, t) 

dt^-i{Sk, t) = J2iHi{Sk, Si)^^i{Si, t) + H_^{Sk, Si)MSi, t)+ 
I 

+ Hi{Sk, Si)iJi{Si, t)) + H_i{Sk, Si)ij^i{Si, t) 

We would now like to describe a corresponding dynamics in the usual space. We recall that 

Sk is both a point in a configuration space as well as a subset of a regular space; namely, a 
set of points x satisfying ui{x) a; k. To make this more precise, we can define Sk as 

Sk^{xe R^\\ui(x) - A;| < e} (10) 

for some small e. In light of the fact that Ui{x) (just like ^2(2;)) can be approximated by 
integers for almost all x, the union of all Sk covers most of our space. The small "gaps" 
that are left as a result of continuity of Ui are not very important. The value of ip{Sk) 
is represented by ipix) for any x & Sk- Therefore, whenever Ui{x) and Ui{x') are both 
approximated by the same integer k, we would like to have 

ipi{x) a ipi{x) ; V-i(^) ~ 'ip-ii^') ; ~ ) ; 'ip-i{x) ^ ^p-,{x ) (11) 

The fields Hi, H^i, Hi and H^i should also be translated into the usual spacetime. By 
inspection, H{Sk, Si) is defined to be a function of two points in the configuration space: 
the one that emits a signal (Si) and the one that receives it (Sk)- Thus, we should identify 
a region in a regular space, Uki, that somehow "encodes" both Si and Sk, and H should be 
nearly constant within that region. That is the main purpose why two fields, Ui and U2 were 
introduced. In particular, we define it to be SkHTi, where 

Ti = {x\u2ix) -l\<e} (12) 

This means that, in order for Hi, H_i, Hi and H_i to be consistently defined, we would like 
to introduce the dynamics of these four fields in such a way that 

Hi{x) ^ Hi{x') ■ H_i{x) ^ H^i{x') ; Hi{x) ^ H0) ; //_,(f) ^ //_,(f' ) (13) 

whenever both ui{x) ~ ui{x') and U2{x) ~ U2{x') hold. Now, once we have established the 
consistency of t/j-s and H-s, the "translation" of the four evolution equations becomes 

dtM^,t)^ [ d^x{Hi{x')^P_i{Si,t) + H_i{^)^Pi{Si,t)+ 
Jti 

+ //,(f )V-i(f', t)) + H_i{^)i^_i{Si, t)) 
dt^_i{x,t)^ I d^x{Hi{x')A{Si,t) + H_i{x')ij_,{Si,t) + 

+ Hi{x')M^: t) + H^i{x')tl;^i{Si, t)) (14) 

dtilJi{x,t)^ [ d''x{Hi{x')ilJi{Si,t) + H_i{x')i/j.i{Si,t) + 
Jti 
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+ Hi{x')^|J_i{x',t) + H_i{x')^|Ji{Sl,t)) (15) 
^t^|J-i{x,t) ^ [ d^x{H,{x')ij.,{Si,t) + H_,{x')t/ji{Si,t)+ 

+ Hi{x')i^i{x\t) + H_i{x')i^^i{Si,t)) 

To sec why it works, consider the first equation. Since there is a one to one correspondence 
between 5*^ H T/ and (-Ui, ^2)1 we know that different Sk H Ti fill most of the space. Thus the 
integral can be approximated as a sum of integrals over each 5'^ n T;. Now, in the interior of 
SkHTi, Hi{x) H{Sk,Si) and i/j-i{x, t) pa ip-iiSi, t). Therefore, if the volume of 5'^ n Ti is 
denoted by V{Sk n TJ), the expression becomes 

dtM^', t)^Yl ^("^^ ^ Ti){Hi{Sk, Si)ip.i{Si, t) + H_,{Sk, Si)iJi{Si, t)+ 

k,l 

+ H,{Sk, Sl),p_^{Sl, t) + H_i{Sk, Si)MSi, t)) 

Now, due to a very large volume of the universe as well as very fast fluctuations of Ui and 
U2, the volumes of different Sk n Ti approximate a common value V. Thus, the expression 
becomes 

dtM^, t)^V J2iHi{Sk, Si)^-i{Si, t) + H^,{Sk, Si)7j;i{Si, t)+ 

k,l 

+ Hi{Sk, Si)^.i{Si, t) + H_i{Sk, Si)MSi, t)) 

which is the same as the original expression we had for configuration space. The other three 
equations work out similarly. 

Therefore, our goal for the sections that follow can be summarized as follows: 

1) Describe dynamics of ip in such a way that 

a) is nearly the same at points with similar values of Ui 

b) ip can be communicated between different points by the integral equations given above 

2) Describe dynamics of H in such a way that 

a) H is nearly the same at x and x' as long as Ui{x) ~ Ui{x') and U2{x) ~ U2{x') 

b) H corresponds to the global Hamiltonian predicted by standard quantum mechanics 
(or standard quantum field theory). 

In the next section we will assume that H is already given, and we will focus on dynamics 
of ■0- Then, in the section after that, we will return to H. 

5. Dynamics of ijj when H is given 

As was mentioned before, in terms of our usual space, the points Sk and Si in configuration 
space are defined as 

Sk^{xe R^l \ui{x) - A;| < e} ; Si ^ {x e \ui{x) - /| < e} (16) 
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for some small e. Our desired dynamical equation represents the communication from Si to 
Sk- That communication happens by means of points x simultaneously satisfying Ui{x) ^ / 
and U2{x) ~ k. The signal emitted by these points is to be received by points x' satisfying 
ui{x') ~ k. Thus, we would like the frequency of emitted signal to be approximately k. This 
can be accomplished by postulating a dynamics that guarantees frequency of emitted signal 
to be U2{x). 

We then encounter a problem: as the signal travels through space, the U2 changes, 
which means that we run the risk of the frequency of that signal changing as well. We will 
simply avoid this issue by coming up with a mechanism of random "sparks" that trigger the 
emission of signals with frequency U2- Since these sparks are localized both in space and in 
time, only the value of U2 at the place of their occurrence is important. After the signal had 
been produced, its subsequent motion is identical to the one in a vacuum, regardless of the 
behavior of U2- Thus, original frequency, 1*2(^0) is preserved. 

The specific mechanism of spark production is not very important. For definiteness, we 
postulate some real valued field x that evolves according to usual wave equation, 

^ - c,V^x = (17) 
and simply postulate that the intensity of a "spark" is given by 

f{x) = l + -tan-\a^(x-Xo)) (18) 

where is a very large number. Thus, /(x) is close to if x < Xo ~ and it is close to 
1 if it is greater than xo + for some small S. Due to the compactness of the universe, 
sometimes different Fourier components of x will produce local resonances. Xo is chosen in 
such a way that it is not likely for x to exceed Xoi thus f{x) looks like a sequence of random 
sparks. At the same time, xo is small enough for these "sparks" not to be too rare. 

Now, we are ready to discuss the emission of the signal. The sources of emission of ipi, 
ijji and ip^i will be denoted as e^^, e^ -^, e^- and e^ -^, respectively. In order to obtain 
our desired equations, we would like the amplitude of these pulses to be determined by the 
couplings of /(x), H-s, and -^-s, where the indexes in the two latter fields correspond to 
the indexes in our desired equations. Since we want e-s to be localized in space, we do not 
include spatial derivatives in their dynamics, thus making sure that they don't propagate. 
Thus, we propose the following equations: 

dle^^{t,x) = -U2{x)e^^{t,x) - Xedoe^^{t,x)+ 

+ fixit, x)){Hi{x)ij_,{t, x) + H_i{x)ij,{t, x)+ 
+ Hi{x)iP^i{t, x) + H_i{x)^i{t, x)) 
dle^_^{t,x) = -U2{x)e^_-,{t,x) - Xedoef_^{t,x)+ 
+ fixit, x))iHiiS)iPiit, x) + //_i(f )V'-i(t, f )+ 
+ Hiix)i;iit, x) + //_,(f)V'-i(i, f) 
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dQe^^{t,x) = -U2{x)e^^{t,x) - Xe^oe^^{t,x)+ (19) 

+ f{x{t,x)){Hi{S)iMt,S) + H^i{S)iP-iit,S)+ 
+ Hi{x)ijj-i{t, x) + H_i{x)ijji{t, x) 

dle^_Xt,x) = -U2{x)e^^^{t,x) - Xedoe^_^(t,x) + 

+ fixit, x)){Hi{x)^^i{t, x) + x)+ 

+ Hi{x)ipi{t, x) + H_i{x)ij_i{t, x) 

In the above equations, the first two terms assure that whenever an excitation occurs, a short 
pulse is sent, with frequency U2. The last two terms tell us that the source of excitation is 
a coupling of /(%), H and ip. The indexes selected in a way that correspond to our desired 
evolution equations of four ip fields, written earlier. As we have pointed out earlier, the above 
equations don't have space derivatives, and thus e-s do not propagate. The communication 
between different points in space is enforced through messenger fields, /i^^, and 
|J,^_^. Their dynamics is given by 

^tV^. - cJVV^. = e^i ; 5tVv-i - cJVVv-i = e^-i (20) 

where represents spacial Laplassian and is some very large constant that corresponds 
to the speed of propagation of /x. Once again, none of the speeds of the signals have anything 
in common with the speed of light; in fact, they are expected to be much larger than that. 
Thus, we have found a mechanism by which the signals are emitted with frequency 1*2(^0)) 
as desired. Since their further propagation is mathematically identical to the one in a 
vacuum, their frequency stays constant. Thus, when that signal passes a point x for which 
U2{x) 7^ U2{xq), its frequency continues to be equal U2{xq). 

Finally, these fields are received by the corresponding "receptors" at the desired loca- 
tions, r^,j, r^_i, and r^_^. As mentioned earlier, the absorption resonance frequency is 
based on ui rather than U2 and the process is defined as follows: 

dfr^^{x, t) = -ui{x)rf^{t, x) - Xrdor^j^{t, x) + n^p^^it, x) 

dtr^_^{x,t) = -ui{x)r^_-,{t,x) - \rdor^_^{t,x) + fi^_^{t,x) 

d^r^.{x,t) = -ui{x)r^.{t,x) - \rdor^.{t,x) + fi^^{t,x) (21) 

dtr^_i{x,t) = -ui{x)r^_Xt,x) - Xrdor^_Xt,x) + ^^_^{t,x) 

The above equations are mathematically identical to the driven harmonic oscillator in r^, 
with sinusoidal driving force fi^. While Ae was assumed to be large enough for the e- 
oscillation to have a short lifetime, we make the opposite assumption about A,.: it is assumed 
to be very small. As a result, whenever ui{x') is close to the frequency of //^, this results 
in a very large resonance in r^. But, the frequency of fi^ is equal to ^2(2;), where x is a 
point of emission of a signal. Thus, whenever U2{x) ^ Ui{x'), x' absorbs signals emitted at 
X. On the other hand, if u{x) and «(£') are different, the fact that both are approximated 
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by integers implies that their difference has to be an integer; that difference is large enough 
for the signal to go virtually unnoticed. 

We now assume that our space is compact (say, a large sphere). Then the emitted signals 
circle our universe multiple times. Noticeably, the "absorption" of fi^ in above dynamics 
is not conservative: it only has an effect on tp, but not on /i^. The emission, on the other 
hand, does have an effect on /x^. As a result, /x^ fields accumulate over time. Most of them 
are "old enough" to have circled our universe multiple times, which means that location of 
their emission is no longer important: they are just as intense around the points of their 
emission as they are billion light years away. 

If their speed is large enough for them to circle the universe within very short time, 
then, on larger time scales, the /i-signals emitted during a physical process have global effect, 
which is independent of location. This is one of the key features that leads to the observed 
non-locality of the configuration space. Noticeably, their accumulation in time leads to the 
increase of ipi, ipi and But, as mentioned earlier, when it comes to the total field 
■0, the increases of these fields pairwise cancel, which results in the expected behavior of ip. 

We know from the theory of oscillations that the squares of the amplitudes of out-of- 
phase oscillations add up. So, since we have established that the amplitude of our oscillations 
is position- independent, we will get the desired dynamics of ipi, ip-i, ipi and ip^i if we associate 
them with squares of amplitudes of oscillation of r^^, r^_-^, r^. and r^_^, respectively. Since 
these are four separate fields, each of these four relations is enforced separately, independent 
of the other three. In other words, contrary to what we are used to, each of the "components" 
of is, itself, a square of the amplitude. 

This can be accomplished by setting up the dynamics of ipi in such a way that whenever 

ijji is less than r^^, it "catches up" fast; but if it is greater than r^^, it decreases slowly. 
Suppose the peak of the oscillation of r^^ is R^^ . Then, as soon as r^^ has reached 7?^^ , the 
ipi quickly reaches the same. But then when r^^ goes back to 0, the ipi "doesn't have time" 
to decrease much, until is back to again. Thus, if the frequency of oscillations is 
high enough, ipi ?s R"^^ at all times. 

On the other hand, if the amplitude of the oscillations has changed from R^-^ to S^-^ < 
R^^ , then after enough time passes, the ■^i field will finally reach S"^^ (even though it might 
take several oscillations for this to occur). Thus, the amplitude of ipi will behave the way it 
is expected on a larger time scales. The same, of course, is true for "0-1) V'i -0-1 . Thus, 
we postulate the dynamics of these four fields to be 

9^01 = a(e'('-V'i-'^i) - 1) ; 9o'0_i = a(e'^'"*-i-'^-^^ - 1) 

9oV, = a(e^(4-^^) - 1) ; dl^l^^, = a{e'^'''^-r^-^ - 1) (22) 

It is easy to see that if a is very small and b is very large, this will lead to the desired result. 

We have just established that the field ■^i is proportional to the square of the amphtude 
of oscillation of r^^. Furthermore r^^ is simply a driven harmonic oscillator with a driving 
force /i^j .Thus, from the theory of driven harmonic oscillator, we know that the amplitude of 
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oscillations of r^^ is proportional to the amplitude of //^^. By inspection our desired equation 
from the previous chapter, 



+ 



+ Hi{x)^l^_^{x,t)) + H_i{x)xl^^{Sut), (23) 

the latter is made up by incoherent signals coming from the collisions of (-0, —i) with [H, r) 
and (-0,-1) with {H,i). 

6. Dynamics for H 

In the previous section wc have assumed that H was only a function of x and not of t. 
Physically, this means that the probability amplitude of transition between any two given 
states is the same, at all times. The only thing that changes with time is their probability 
amplitudes. However, while this agrees with the standard view, we are not happy with 
that view. After all, we don't like the fact that transition probabilities between different 
configurations are non-local. Therefore, we would like to come up with a local mechanism 
by which H attains this global value. We have shown a jist of it in section 3, when we have 
discussed the difference between h and H. 

The good news is that the sought-after mechanism should be based on the configuration 
spaces alone; that is, the configuration of the particles. As was explained in previous sections, 
particles are stationary. Therefore, the task of obtaining global equilibrium for H is much 
easier than the one for ip: H can evolve as slowly as we wish, since it doesn't have to "keep 
up" with any changes. The challenge, however, is that, at any given point with Ui{x) ~ k 
and U2{x) ^ /, we would like H{x) to look only at the states |si > and \s2 >. 

In order to be able to "filter out" the rest of the information, we will have to introduce 
harmonic processes with frequencies ui and U2 at every point. Formally, this can be done 
by introducing two fields cti and a2- In order not to "confuse" the oscillatory processes 
associated with different points, we should make sure these fields do not propagate. Thus, 
their dynamics does not have space derivatives. We postulate their equations of motion to 
be 

d^ai{x,t) = —ul{x)ai{x,t) ; d'^a2{x,t) = —ul{x)a2{x,t) (24) 

The above equations, however, allow for solution to alter in space both its phase and ampli- 
tude: 

ai{x, t) — Ai{x) cos{t — ti{x)) ; a2{x, t) — A2{x) cos{t — t2{x)). (25) 

Since we do not want the above variations in amplitude, we introduce another two fields, Pi 
and ^2, which are the normalizations of ai and 0:2 respectively: 

y/ul{x)a'l{x,t) + (dtai)^(x,t) y/ul{x)al{x,t) + {dta2Y{x,t) 

Of course, the phases ti{x) and t2{x) can not be dealt with this way; but as we shall see, 
they will eventually become unimportant. 
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Now, in order to introduce the Hamiltonian, we have to first come up with "local" way 
of defining it point-wise and then "integrating" it. Since particles have zero size, they don't 
"touch" each other. Thus, we do not have literal vortexes. Instead, we need to define a 
"density field" associated with each particle that forms a cloud around it. While the particle 
is still a point object, that density field is extended in space. Thus, the overlap of the density 
fields of, say, electron, neutrino, and W-boson, forms euW vortex. Of course, this calls for 
rewriting quantum field theory in a more continuous way. This is far beyond the scope of 
this paper. 

As far as this paper is concerned, we will assume that we do have a version of quantum 
field theory that allows particles to have size, and our only goal is reinterpreting phase space 
and introducing determinism. However, there is one important thing to address: if a particle 
has size, it will undoubtfuUy overlap with different S-s. But we want that particle to be 
detected only in the Sui{xi)) and 5'„2(xo)) where xq is a location of a "center" of the particle. 
We will do that by introducing the emission and reception mechanism of p-field. That field 
is being "emitted" at xq and then "received" at x that belongs to the same Su- 

Let us now be a little bit more explicit. Suppose we have A types of particles. For 
example, particle number 1 is electron, particle number 2 is proton, etc. We will denote the 
density field corresponding to particle number A hy pa- The "source" of that p field is a 
"point charge" qa, and its "messenger" is pp^. The latter undergoes oscillations, which are 
guided by the behavior of /3 at the center: 



where k is the numbering of different particles of the same type A (they are assumed to 
be distinguishable), and XAk is their location in space (there is no ^-dependence of XAk 
because, as explained earlier, we have assumed they are stationary), in the above expression 
time delay is introduced for the purposes of locality. The inverse dependence on distance 
assures a very small size of density distribution, given the appropriate value of qa- Since any 
resemblance with electrostatics is merely the issue of convenience, any differences between 
the two theories are not very important. 

Now, at any given point, we would like to define a local process. For example, suppose 
near a given point there is an electron, a neutrino and W boson. We would like to say, 
for example, that what contributed to Hamiltonian is the electron emitting W boson, and 
turning into neutrino. Thus, we would like to view electron as part of Sui(x), while viewing 
both neutrino and W boson as parts of Su2{x)- This can be accomplished if electron, neutrino 
and Vt^-boson are placed at points xi, X2 and X3, all of which are close to xq and which 
satisfy Ui{xi) — Ui{xo), 1*2(^2) = 1*2(^0) and 1*2(^3) = ■^2(^0)- Thus, a point Xq needs to 
have receptors for both ui and U2- The same applies to every other point. Thus, we define 
two separate receptors: 




(27) 




/^pa(^> t) - ul{x)rp2{x, t) - Xr^ 



Pp^{x, t) - ul{x)rpi{x, t) - \r, 



PA2 



dt 



dr 



PA2 



PAl 



(29) 



(28) 
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As before, the receptor fields oscillate; but we do not want our p to do that. We use the 
same tactic as in previous section, and " convert" receptor fields into non-oscillating densities 
via the following expressions: 

^ = a.le^^C-?..--) - 1) ; ^ = a,(e^^(^^..--) - 1) (30) 

We are now ready to move on to define Hamiltonian. As stated before, we are dealing with 
modified version of quantum field theory (which we never introduced but we pretend that 
we did) that is based on densities of the particles as opposed to particles themselves. Thus, 
if the above local neighborhood was all there is in the universe we would simply say that 
H is a function of p-s. But since that is not the case, we would have to "integrate". The 
way to do that locally is to follow previous procedure and introduce emitter (e//), messenger 
field {ph) and receiver field {th)- Thus, only the emitter field is the local function of above 
described densities. 

However, given that our emitters are placed fixed distances away from each other, if the 
emission was constant that would lead to interference. Since we do not want that, we would 
like to replace a continuous process with discrete emissions. When emissions do occur, they 
have prescribed frequencies; but, at the same time, they are pulse-like. If we make sure that 
they are random in time, the total phase shift between emissions at any two given points 
will be random, despite the fact that the distance is known. This will assure us that, on 
average, the squares of amplitudes add up, regardless of the distance. 

In order to come up with deterministic mechanism of producing "random pulses" we 
need a "random generator". For that purpose, we introduce another field, x(a;,t), that 
evolves according to usual wave equation, 

dh - clV\ = 0. (31) 

As with all other fields, is much larger than the speed of hght. In fact it is so large that 
within small time interval, the field x circles our universe multiple times (as stated before, 
we assume that our universe is compact). Since the exact shape of the universe is not known, 
there are might be some resonances of x random places. We now introduce a differentiable 
approximation to step function, 

/(^) = ^ + ton-i(^-^o) (32) 

If we choose in such a way that x(t, x) exceeds 9q only at very unlikely resonances, it will 
mean that /(x(t, if)) is 1 at these rare instances and everywhere else. In other words, the 
latter will look like a set of randomly distributed pulses, which is what we need. Thus, we 
define our emission equations for ifi, and H_i as follows: 

dltH^x^t) = -uj'^{ui{x),U2{x))eHi{x,t) - XedQeHi{x,t)+ 

+ f(x(x, t))hi(pn{x, t), pbi(x, t);pi2{x, t), pb2(x, t)) 
^e.H_^{x,t) = -u'^{ui{x),U2{x))eH_Ax,t) - XedoeH_A^,t)+ 
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+ fixi^, t))h_i{pii{x, t), pbi{x, ty,pi2{x, t), pb2{x, t)) (33) 
dteHi{x,t) = -uj'^{ui{x),U2{x))eHiix,t) - Ae9oeij.(f , i)+ 
+ f(x(x, t))hi(pii(x, t), Pbi(x, t);pi2(x, t), Pb2(x, t)) 
dfeH_i{x,t) = -uj'^{ui{x),U2{x))eH_i{x,t) - XedoeH_iix,t)+ 

+ f(x(x, t))h-i(pn(x, t), pbi(x, t);pi2(x, t), pb2{x, t)) 

In the above expression it doesn't matter that and Cij. are simultaneous both in space 
and time. After all, from the previous chapter we remember that we did not directly add 
ipi, ip-ii ipi and ip^i until the very end, at which point we no longer had to worry about 
their oscillations. Similarly, we do not add and either until that point. Thus, the 
only thing we need to worry about is the interference of these fields with themselves, and not 
with each other. 

Since f{x{x, t)) is either or 1, the amplitude of the pulses are determined by h. Thus, 
losely speaking, after the propagation, h'^ will globally " add" to give us H. The mechanism 
of this propagation is similar to the one of ip. We introduce the "messenger fields" /iHr 
II Hi which are subject to the following dynamics: 



o2 2 V72 I \ 



^'^^^H. - cl^^ V Vh. + ^ = (34) 



o2 2 vt2 I \ '-^H'tt-i 

OtHH.i - c^^_ V p,H_, + Kh^^ = eH_, 



„2 V72. , X d^lH^ 

dt 

Then, by the same reasoning as with ip, we claim that if we introduce "receptor fields" 
they will be able to "listen" to everything emitted in the entire universe (since the speed 
of propagation of signal is very large, the time delays are negligible), which is the ultimate 
source of non-locality. Similarly to previous section, we define them based on driven harmonic 
oscillators, as follows: 

dlrHi{x,t) = -u'^{ui{x),U2{x))rHi{x,t) - Xr„^dtrHiix,t) + iJ,Hiix,t) 

d^rH_,ix,t) = -uj'^{ui{x),U2ix))rH_^ix,t) - Ar^_^(9tr^_, (f , t) + fiH.iix,t) (35) 

d^rHi{x,t) = -u'^{ui{x),U2{x))rHi{x,t) - A^^, ^tr^, (x, t) + fj,Hi{x,t) 
dtrH_i{x,t) = -uj'^{ui{x),U2{x))rH_i{x,t) - \rH_ArH.A^,t) + fJ^H_i{x,t) 

Finally, we would like to get rid of unwanted oscillations of by defining H to be an 
amplitude of the above oscillations. By the same arguments wc used in the previous section, 
the way to do that without violating locality in time is through the following dynamics: 

d^H^ = a(e''(''^i-^^) - 1) ; d^H_, = aie'^'"^-^-"-'^ - 1) (36) 

d^Hi = aie'^'h-''^^ - 1) ; d^H_, = a{e''^''"-r"-^^ - 1) 
where b is very large. 
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7. Beables 



7.1 Review of Diirr et el 

Up until this point we have shown how to model quantum field theory amplitudes in terms 
of functions only over R^. In this chapter we will do the same with theory of measurement. 
We will follow similar philosophy in a sense that we will start out with global theory and 
then "convert" it into local one. In principle this can be done with any pilot wave model. 
However, in light of the fact that our configuration space is discrete, it is easiest to use the 
model of stochastic jumps due to Diirr et el (see ^). Other models, on the other hand, would 
likely have to be artificially discretized before they can be "plugged into" our framework. 
For that reason, in this paper we will focus exclusively on Diirr et el, while leaving the rest 
of the Pilot Wave models for future work. In this subsection, we will briefiy outline the 
global picture proposed by Diirr et el, and then in the next subsection we will move on to 
"encoding" it into M^. 

The key concept of Pilot Wave models is that a particle and a wave are completely 
separate substances. There is no such thing as " collapse of wave function" . Instead, a wave 
evolves according to Schrodinger's equation (or some other rules of quantum mechanics) at 
all times. At the same time, a particle is local at all times, and it is being guided by a 
wave according to guidance equation. In case of non-relativistic quantum mechanics, that 
equation is 

= ±^ ImlmP, (37) 
at m 

but that equation might become increasingly more complicated for more general cases, specif- 
ically involving quantum fields. In all cases, however, the guidence equation has to be de- 
signed in such a way that if we don't have enough information about the particle, then the 
probability of finding it will end up being {ipl"^- 

Let us now translate the above into the language of multiparticle states. The notion of 
wave function generahzes to an assignment of probability amplitudes to all possible states. 
Schrodinger's equation generalizes to Feynman rules. And, finally, the notion of localized 
particle generalizes to one specific state. Thus, the claim of any Pilot Wave model is that, 
at any given time, we assign the probability amplitudes to all states, and, at the same 
time, we have one specific state that evolves in time, guided by these amplitudes. 

One of the main challenges of this, however, is that the evolution of state involves 
creation and annihilation of particles, which is not a continuous process, while the above 
guidance equation is continuous. Diirr et el decided to address this issue by replacing the 
velocity v with a probability of a jump, a{e, e') from the state e to the state e'. This, 
however, violates determinism because, while we know the probability, the actual timing of 
each jump is random. Thus, in order to be as close as possible to determinism, they have 
claimed that between the jumps our state evolves according to differentiable (and, therefore, 
deterministic) equation. That process is simply being "interrupted" by the discrete jumps. 

Now, in order to calculate the probability of these jumps, we have to compute the time 
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derivative of desired probability of each state, based on quantum field theory. After that, a 
suitable Pilot Wave model is introduced in order to accommodate that desired probability. 
Now, suppose the evolving state is >, and the probability of a state \e > is ■ip{e) =< el'if > 
Then the total current flowing into that state is given by 

|(,.» = 2«e(.-^) (38) 
Now, from evolution equation we know that 

^ - i < e\H\^ >=iJ2< e\H\e' >< e'|^ > . (39) 

e' 

Combining this with 

i/j* > (40) 

we get 

-^(^*^) = 2ImJ2i< >< e\H\e' >< e'\^ >) (41) 

e' 

This means that we would like a current from |e' > to |e > to be given by 

j(e', e) = 2 Im{< >< e\H\e' >< e'|* >) (42) 

Now, that current is equal to the probability of a jump, multiplied by the probability density 
of e' (namely, ip*ip =< e'|^ >< ^^je' >). Thus, in order to obtain the probability of the 
jump, we have to divide the above expression by the latter. Thus, our first guess is 



^(^'^ =• 2 Im ' . (43) 



Naively, this would give us a negative probability of jumps in direction opposite to the 
current. This, of course, wouldn't make sense. So, instead, we claim that the probability of 
the latter is 0. We can do that aby introducing a function which is equal to x when x > 
and when x < 0. This gives us 



cr(e , e) = 2 Im ' — ' — — —r- — — (44) 



< * 


\c >< c\ 


\H\ 


\c' >< c'\ 


> 


< e'l 


^ >< * 


\e' > 



(< *|e >< e\H\ 


|e' >< e'l 


|*>)t 


< e'l 


^ >< * 


|e' > 



As we previously said, however, according to Diirr et el, the above is just part of a dynamics. 
In particular, according to that model most of the time the system evolves according to a 
differentiable, deterministic equation. It is merely being interrupted by stochastic jumps with 
the above probability. However, Diirr et el also notices that continuum process, itself, might 
be simply a limit to the discrete. In particular, they quote Bell's suggestion of representing 
fermions in terms of the numbers of particles at every discrete lattice points. 

The ultimate reason why Diirr et el tends to favor the option that does include continu- 
ous processes over the time intervals is, obviously, because it is more deterministic, although 
still not completely so. In our work, however, we will show how the jumps themselves can 
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be modeled as deterministic processes. In the language of introduction, we are building a 
"computer" that "simulates" anything and everything we might want. We know that a com- 
puter is ran by different iable, deterministic processes. Thus, as long as we will successfully 
simulate something, it will be continuous and deterministic, by default. In other words, the 
jumps are no longer "worse" than the continuous process. 

On the other hand, in light of the fact that our configuration space is discrete, model- 
ing continuous process is very problematic; attempting to do so would probably require to 
discretize it first! For that reason, we will assume a Pilot Wave model that is being ran by 
above-described jumps alone, without any continuous evolution. That, of course, is possible. 
After all, we don't have access to very small scales so whatever appears continuous to our 
eyes might, in reality, be discrete. 

7.2 Converting Diirr et el into a local theory 

We are now ready to "convert" the non-local model of Diirr into our "local" theory. In order 
to do that, we will introduce one single beable particle, whose position is x = b(t) G M^. 
We will also introduce beable field S: x R — )■ M. Our goal is to find a dynamics that 
would guarantee us that B{x,t) is non-zero only when ui{x) = ui{b{t)) and, in the latter 
case, B{x,t) is roughly independent of the choice of x. 

If the above is done successfully, it guarantees us that the position of single particle 
in R'^ (which we will call b-particle) "encodes" a configuration of multiple particles that is 
realized at a given time. In particular, the b-particle is like a "lamp" that "illuminates" the 
particles that have approximately the same value of ui, and doesn't illuminate anything else. 
Then, when our "lamp" moves, the values of ui the places it occupies change. Thus, the 
set of particles it "illuminates" changes as well, which leads to an illusion that they move 
around. 

We would like to be able to find a local theory that imitates non-local behavior described 
in Diirr et el (^). In that paper, they claimed that within small interval dt, the probability 
of a "jump" between the state |e > and |e' > is a{e, e')dt, where 



< ip\e >< e\H\ 


\e' >< e'\ip > 




\ <ip\e> \ 


2 



(45) 



In order to mimic that result, we have to claim that, over some short period of time, the 
probability that b-particle would travel from x to x' is proportional to a{Su(x), Su{x')), where 
is a quantum state that can be "read off" from u ^ uq. But then the question is: what 
is a local mechanism for b-particle to "scan" the values of u{x') at various x' and "decide" 
where to move? 

Thankfully, from previous sections, we already have a clue of where to look at: we 
have two fields Ui and U2- We recall that, as far as quantum mechanical amplitudes were 
concerned, every point x can serve as a "transition" from the state Suj{s) to the state Su2{x)- 
We will now introduce similar concept for beables. In other words, the point x' to where 
the b-particle transitions to, has to satisfy ui{x') ~ U2{x). The reason ui{x') appears to be 
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random is that wc only know Ui{x) and not U2{x). After all, only the former can be "read 
off' from the configuration of the "visible" particles. 

Again, the above is non-local, so we need to come up with some local mechanism for 
that. Furthermore, we also have to come up with a way in which a particle "decides" to 
which of the x' e S^^^g) to transition to, which trajectory to take, etc. The mechanism that 
we propose is that a particle sends off signals that cause a " potential" to fall throughout the 
entire Su2{x)- After that, it is attracted to that same potential that it caused. Its specific 
trajectory is determined based on its initial location. 

Before we proceed with this, we have to make sure that we enforce the probability a{e, e') 
defined above. We will first find a local way of defining the desired probability, cr(a;), and 
then afterwards propose a deterministic mechanism for that desired probability to be the 
actual one. In order to come up with desired probability, the information about %Ij{x') (or, 
in other words, ip{Su2{x))) has to be available at x. 

This can be done by furnishing point x with an "antenna" tuned to U2 rather than ui, 
and thus introducing another receptor fields, r^i and if)' that are being "read" from that 
antenna. Since the mechanism of obtaining r^i{x) and ^^'{x) is the same as the one for r^{x') 
and ip{x'), respectively, the former should approximate the latter. The good news is that 
the former represents fields at x, and not at x'. This allows us to come up with a local 
expression for a{x,x'). 

Now, the mechanism of getting r^r and ip' is a carbon copy of the equations for r,p and 
ijj, where ijj is replaced with ip', r^ is replaced with r^r and Ui is replaced with U2- However, 
we do not replace with e^/; and we do not replace /x^ with /x^/ either. The reason for this 
is that we would like r^i{x) and r^{x') to be "listening" to the same thing, namely /i^ (and 
not /i^i) in order to be equal. 

As before, we need to have four sets of equations corresponding to indexes 1, —1, i and 
—i. Thus, the set of equations of r^/ is 

d^r^'^ {x, t) = -U2{x)r^'^ {t, x) - A^^or^'^ {t, x) + /x^, {t, x) 

d^r^'_^ (f , t) = -U2{x)r^'_^ {t, x) - Xrdor^>_^ {t, x) + {t, x) 

dtr^',{x, t) = -U2{x)r^',{t, x) - \rdor^>.{t, x) + //^,(t, x) (46) 
d'^r^i^^{x,t) = -U2{x)r^>_^{t,x) - \rdor^'_^{t,x) + fJ.^_i{t,x) 
Our way of obtaining ip' from r^i is the same as the way we obtained from ip, which is 

c^V; = a(e'"*r*» - 1) ; a„V^ = a(/"5-.-*- - 1) (47) 
Now, as we have just mentioned, we interpret the values of ip and ip' as 

^Pix) =< ij\Su,{x) > ; i^'ix) =< ij\Su2{x) > (48) 
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Furthermore, as we explained in previous sections, the Hamiltonian is encoded in H{x) via 



< Su,{x)\H\SuM >= H{x,t) 



(49) 



Substituting these into the equation proposed by Diirr et el, we obtain 



(^{t; Su2{x)) = cr{x, t) 



(50) 



where 



a{x, t) 



Im{ip{x, t)ip'{x, t)H{x,t))^ 



(51) 



The above gives us the desired probability of the transition of the b-particle. Now we would 
like to discuss the mechanism of that transition as well as how that desired probability is 
enforced. As we said before, we would like our b-particlc to emit a signal that will cause the 
potential to fall throughout Su2{x)- Thus, we would like to find a deterministic mechanism 
of emission of a signal in such a way that it appears to occur with probability a{x,t). We 
do that by introducing a-particles whose density is proportional to a{x,t) (while there is 
only one b-particle there are many a-particlcs), and the signal is emitted if and only if our 
b-particle collides with one of the a-particles. 

Since the particles have zero size, they have zero probability of collision. Thus, we have 
to define a field p that surrounds each particle, and by "collision" we really mean the overlap 
of these fields (intuitively, this is equivalent to saying that particles have size, and p is the 
density of matter distribution that makes up these particles). We will denote these fields, 
corresponding to particles a and b, as Pa and pb, respectively. It doesn't matter what specific 
expression we choose for them as long as it meets the above property. For definiteness, we 
will postulate 



Here, we put a summation sign for a-particles but not for b, because we assume that there is 
only one b particle, which corresponds to a single state defined by a beable. Now, in order to 
enforce the desired density of a-particles, we come up with a dynamics that assures us that 
their velocity is inversely proportional to the desired density. Thus, the larger the desired 
density is, the slower they move, and thus they spend more time in that region. In order for 
the direction of their velocity to be specified as well, we will use acceleration to attain the 
desired velocity. The direction of that acceleration is parallel to direction of velocity at any 
given time: 



In other words, from the space-alone perspective, a-particles always move along geodesies, 
with varying speeds. But in light of the compactness of the universe, these particles come 
back moving in a very different direction after having circled the latter. This provides 
mechanism of randomness of their velocity directions. 





(52) 




(53) 
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We will now have to introduce a mechanism of generating the desired potential when 
the overlap between pa and pf, occurs. As always, we will do that through emission field ey, 
messenger field /xy, reception field ry and finally a potential V . This is a bit tricky though. 
Suppose b-particle, after collision of a-particle, was "sent" to some region Su2{x)- But on 
its way there it was crossing a region Su^, and there it underwent collision with another 
a-particlc that sent it to S'^^. In order to avoid this situation, wc have to make sure that the 
"source" of ey is large only if the velocity of a particle is small (thus, the particle can only be 
"sent" to a new destination after it has "arrived" to the previous one, and, therefore, slowed 
down). Since the interaction occurs through overlap of p-fields rather than actual collision 
of particles, the "velocity" relevent to this interaction is the "current" of pb, defined as 



The reason we have only introduced jf, and not ja is that our goal is to avoid the interaction 
when b-particle is in transition. The latter is characterized by jb being large. The values 

— * 

of ja are independent of that and, therefore, should be irrelevant. In the above expression 
we used Cp instead of coming up with a separate constant Cj, and that is for a good reason. 
If the values of Cp and Cj were different, then the two fields would have different ranges of 

— * 

interaction. Now, if the range of was smaller than the range of pi,, then some distance 
away from b-particle it would appear that the latter is stationary. Thus, in that region a-b 
interaction will occur "by mistake". Since we don't want that, we make sure that Cj^ = Cp^. 

Apart from wanting emission to be small whenever is large, we also would like it to 
be small when the magnitude of a "potential", |V"|, is large, as well. Suppose a particle 
makes a transition from xi e Su^ to X2 & and then, later, it makes another transition 
from X3 G to X4 G (the reason we have four x-s is that a particle can move while in 
Su2, as will be discussed later). Now, the fact that the transition from Xi to X2 occurred, 
indicates that there was a fall in potential throughout Su2- Now, suppose that potential 
"didn't have time" to return to after the second transition. Then, the fact that we still 
have potential well throughout Su2 might lead to a particle transitioning to x^ G Su2- We 
do not want that. For that reason, we would like b particle to "wait" until the potential 
becomes approximately before making another transition. Thus, we would like |y| to have 
damping effect on emission, on top of the dumping effect of It is important to indicate, 
though, that even though |V"| is approximately throughout Su2, it is not exactly yet; thus, 
particle still spends most of the time in that region. 

We are now ready to describe the emission field ey produced by a-b interaction. As be- 
fore, since emission field is non-propagating, its equation does not involve any space deriva- 
tives. But, at the same time, in order to avoid 5-function singularities it is defined as a 
function both of space and time. Its equation of motion is the one of dumped harmonic 
oscillator, whose source is proportional to PaPb- In order for the source to be small whenever 
either jb or \V\ is large, we divide it by 1 -|- Nbjb + Nv\V\, where Nb and Ny are some large 
constants. This gives us 




(54) 



d^ev{t,x) 



U2{t, x)ev{t, x) - Xeydtev{t, x) + 




(55) 
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In order to obtain the desired fall in potential throughout Su2{x)^ we would like to introduce 
a messenger signal, iiy that is being emitted by ey. Its equation of motion is given by 

dlixv{t, x) - cl^V^nvit, x) + Xj^yi^vit, x) = iiv{t, x) (56) 

Finally, the desired value of V is produced through the "reception" of that signal. For the 
same reasons as explained in previous sections, the reception is a two step process that 
involves the oscillating reception field ry and then its conversion into non-oscillating V: 

d^rv{t,x) = -ui{t,x)rv{t,x) - \rydtrv{t,x) + iJ,v{t,x) (57) 

d^V{t, x) = kv{e^^^''^-^^ - 1) (58) 

In the above expressions, A^^^ and ky have to be sufficiently large, so that after a transition 
of the b particle is completed, the potential will disappear soon enough; after all, we don't 
want our particle to "wait" for the next transition for too long, since that would distort the 
stochastic nature of the desired jump. At the same time, however, these coefficients should 
not be too large either. After all, if the potential "dies out" too fast, the b-particle might 
not complete a transition. Since its path is continuous, it will be "stuck" at some point x" 
for which u{x") — u^. Given that the ?i-fickl fluctuates quite fast, ^3 will have nothing to do 
with either ui or U2, which would lead to a transition to unwanted quantum state. 

Thankfully, it is quite easy to make sure that the potential dies out neither too slowly 
nor too fast. As we have explained in previous sections, the u field has to be designed in 
such a way that small enough neighborhood of every point intersects every single u^. If such 
is the case, it will take very short time for a particle to make a transition. Therefore, even if 
potential dies out very fast, it might still not be too fast. Thus, the answer to our dilemma 
is for the above mentioned constants to be very large but, at the same time, be smaller than 
some power of the inverse of the scale of fluctuations of u (the latter is much larger than any 
large number we would normally think of). 

There is one more thing to take care of. Since the regions and do not "touch" 
each other, the b particle will not be able to find out a direction in which it has to travel 
based on VV^. We will enforce the desired transition by making sure that the velocity of 
the b particle is large as long as \V\ is small. Thus, a particle will keep moving, until it 
"accidentally" reaches a region of large \V\, where it will "stop". Thus, we would like the 
velocity of the b particle to be 1/(1 + MyV), where My is some very large number. In order 
for the direction of the velocity to be well defined, we will use our earlier trick and impose 
a condition on acceleration; namely, that a particle accelerates in a direction parallel to its 
velocity until its speed reaches a desired value: 



dt^~ ''di\l + NvV 



(- 



dh 
di 



(59) 



The above completes the description of dynamics of b-particle. Thus, we have found a way 
in which b-particle "jumps" into the desired Su^. But it only occupies one point in that set. 
We would like to propose a mechanism by which from that one point it "illuminates" the 
entire Su,^, thus making all of the "electrons" (and other stationary particles) that happen 
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to reside in visible. In previous versions of this paper, we have done this by introducing 
a new field, B, which is produced by a 6 particle. However, right now I realized that the 
potential V has essentially the same properties; thus, based on Ocam's razor, it is better to 
simply use V. But, in order for the notation to be consistent with previous versions of the 
paper, we will still use B and simply define it as B{x) — V{x), This does not agree with 
more complicated definition of B in previous versions, but it serves the same purpose. 

Lets finally explicitly write the way B field "illuminates" our states. As we explained in 
the beginning, the particles we are ultimately interested in (electrons, protons, etc) are all 
stationary. They fall into some regions Suf. , and they are either visible or invisible, based on 
the value of ijj in that region. In light of the fact that we have introduced Pilot Wave model, 
we have to replace ip with B as the ultimate definition of reality. Thus, our particles are 
either visible or invisible based on value of B. Now, as we mentioned before, each particle 
has two short-acting p fields: for a particle of type k they are pki and pk2- They represent 
a "density" of the version of that particle with a size. These fields are independent of how 
much (or how little) a particle is being "illuminated" by B. Now, in order to take that 
illumination into account, we will postulate that the "observed" density of the particles is 
Pk, and is defined to be 

Pk{x, t) = pki{x)B{x, t) (60) 

On the right hand side of the above expression I have chosen p^i instead of pk2 because for 
a particle located at x ^ Sid Tj, pi is felt withing the S'j-neighborhood of x while p2 is felt 
within Tj-neighborhood of x. Thus, in light of the fact that the intensity of B field is a 
function of S-s as opposed to T-s, pi is more appropriate choice. Intuitively, p^i is a sum of 
the finite-width versions of 5-functions, and the S is a coefficient in front of the sum that 
"picks out" which of the terms we "see" and which we do not. In principle, a reader that 
prefers the situation where particles have size is free to replace p^i in the above expression 
with iteral (5-functions. 

8. Conclusion 

In this paper we have presented a way to "encode" quantum field theory amplitudes into 
three dimensional space. Then, later, we used similar techniques to "convert" Pilot Wave 
model due to Diirr et el from configuration space to ordinary as well. Apart from getting 
rid of non-local concept of configuration space, we found a way to make sure that the "jumps" 
postulated by these authors are continuous, and, therefore, deterministic. 

One obvious short coming of this work is that the approach is very much "forced". We 
basically know the kind of answers we want, and we design a "machine" that simulates these 
answers, while reminding ourselves all along that a classically- working " machine" is local by 
default. In some sense this is cheating since any physical concept we "don't like" can be 
explained away by the "machine" that simulates it. Besides, this work implies a fractal-type 
structure: the quantum arises out of classical while it is still true that on larger scales the 
classical arises out of quantum. This might lead to a question: why the coincidence? 

On the other hand, however, one might argue that whenever we conceptualize anything 
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we are doing that based on the concept we are used to. In fact, when we think of space, we 
imagine that we have eyes to see it; or when we think of time, we imagine some changes in 
our brain that make us feel time. Thus, the question "why did fractals arise" can be asked in 
that level, too. It is also important to note that if a cave man was to see a computer, without 
knowing what it is, he, too, might think of some "quantum" concepts similar to configuration 
space; and, if someone were to tell him about the complicated circuits involved, he would 
also say that this explanation is "forced". But that doesn't make him right. 

There are also some unresolved issues on a more technical level. In particular, this work 
was using an "imaginary" version of quantum field theory, where particles have size, and the 
"vortexes" are replaced by the functions of the "matter density". No attempt to formulate 
such theory was made. Instead, it was simply assumed that the configuration space version 
of that theory exists, and we were focused on converting that imaginary theory into the 
language of ordinary space. For the future research, it is important to fill in that gap. 



A. Appendix: Driven harmonic oscillator in classical dynamics 

Throughout this paper, the concept of resonances was extensively used that is simply bor- 
rowed from the dynamics of classical harmonic oscillator. Therefore, for the convenience of 
the reader, in this section we will review the standard treatment of driven harmonic oscillaor. 
The same can also be read in a number of classical dynamics textbooks, including [8]. 

Consider a harmonic oscillator with internal frequency Uq, driven by sinusoidal external 
force of frequency u. The equation of motion of this oscillator is given by 

3^ oc dx 
m-— = sin {ijj{t — to)) — mujQX — A—. (61) 

(JiL (JiV 

In order to solve the above equation it is easier to first solve 

m^l = e^-(*-*o) _ niu^x - (62) 
dt^ ^ dt ^ ^ 

and then take an appropriate linear combination of its solutions. This equation can be 
re-expressed as 

= Ae'^' - m^ulx - A^ (63) 

dt^ ^ dt ^ ^ 

where A = e"*"^*". We will try a solution of the form 

X = Be'""' (64) 
By substituting it into above equation, it becomes 

- mu'^B = A- mulB - iuXB (65) 

which gives us 

m\LO^ — UJ^) + lUJA 
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and, therefore, 

By substituting A = e"""^*"^, we obtain 

m(a;Q — uj^) + luiX 

We now return to the original equation, 

d'^ X dx 
m-— y = sin {uj{t — ^o)) — mujQX — A— (69) 

The solution to this equation is given by linear combination, 
where Xi and ^2 satisfy the following equations: 



m- 



it^ oV-A^ (71) 

and 

= e--(*-*°) - ma;o'a;2 - (72) 

Thus, can be read off from the solution we have just found, while X2 can be obtained by 
replacing uj with —uj. Their linear combination gives 

2i \m{u!Q — ui^) + iuX m{u!Q — u"^) + iuXJ 
which, after some algebra becomes 



m{LjQ — cu"^) sin {u!{t — to)) — coX cos {u!{t — to)) , , 

m(cu2 - cu2)2 + UJ^X^ ^ ' 



By using the identities 

_i LuX \ m.{u!Q — uj^) 



cos tan 



uX \ LOX 



• / -1 '^■^ \ /rrrN 

sm tan — r— -^^ r- = — (75) 

V W?{ujI-UJ^)J ./m2r(w2 ^;2A2 , ^2X2 



the expression for x becomes 



cos [tan ^ 2\ ) sin iuit — tn)) — sin [tan ^ ,%'^2\ ) cos iuit — to)) 
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By using the expression of the sine of the difference, this evaluates to 

^ ^ sin (ujt -to- h)) 

where 

1 _i ^'-^ 

ti = —tan „, 57- (78) 

uj m^[ujQ — uj^) 

The above solution, while correct, is not complete. After all, we can easily choose initial 
conditions on position and velocity that do not meet the above equation. This problem is 
fixed by adding a complimentary solution corresponding to a free oscillator. 

Prom linearity, it follows that if x satisfies the equation of driven harmonic oscillator, so does 
x + Xc- The solution that we have obtained earlier is called particular solution, and denoted 
by Xp. Thus, a general solution is of the form 

X — Xc + Xp (80) 

To obtain Xc, we again look for the solution of the form 

Xc = be"-* (81) 

where a can be any complex number (which is why we put e""* instead of e*"*). Substituting 
this into the above differential equation, we get 

ma^ + Xa + mul = (82) 

which gives us 

If we assume A < 4/c, this becomes 

-A±V4m2a;g-A^ 

2^ ■ ^^^^ 

Since we would like Xc to be real, we express it as a linear combination of the two solutions: 

V . Y g. 



which easily evaluates to 



Xc = Xcoe 2 



-|(.-M,„,(xSEZ(t_t„)). (86) 
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Thus, the complete solution is 

sin (uj(t — tn — ti)) x,. , N / \/Am'^Un — , ,\ 
X = ^ ^ ^ ^ a;coe-2 cos ^ (f -to)), 87 

where 

1 _i tjA 

ti = —tan „, ^ ^ (88) 

a; m^lujQ — u^) 

However, due to e^^(t — tco) factor, the complimentary solution dies out in time. Thus, after 
enough time passes, only particular solution, 

X, = (^(^-^0-^^)) (89) 

survives. The violation of time reversal symmetry is due to the fact that the velocity with a 
minus sign (as opposed to plus sign) enters into the acceleration. Notably, if A is very small, 
the amplitude corresponding io u = is very large. This phenomenon is called resonance, 
and ujq is called resonance frequency. 
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